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Abstract 

Superconducting cosmic strings naturally emit highly boosted charge carriers from cusps. This 
occurs when a cosmic string or a loop moves through a magnetic field and develops an electric 
current. The charge carriers and the products of their decay including protons, photons and 
neutrinos, are emitted as a narrow jets with opening angle 9 ~ l/7o where 7 C is the Lorentz factor 
of the cusp. The excitation of electric currents in strings occurs mostly in clusters of galaxies, 
which are characterized by magnetic fields B ~ 10~ 6 G and a filling factor Jb ~ 10 -3 . 

Two string parameters determine the emission of the particles: the symmetry breaking scale n, 
which for successful applications should be of order 10 9 -10 12 GeV, and the dimensionless parameter 
i c , which determines the maximum induced current as J max = ic e V an d the energy of emitted charge 
carriers as e x ~ icIcV^ where e is the electric charge of a particle. For the parameters rj and B 
mentioned above, the Lorentz factor reaches "y c ~ 10 12 and the maximum particle energy can be 
as high as 7 c ry ~ 10 22 GeV. The diffuse fluxes of UHE neutrinos are close to the cascade upper 
limit, and can be detected by future neutrino observatories. The signatures of this model are: very 
high energies of neutrinos, in excess of 10 20 eV, correlation of neutrinos with clusters of galaxies, 
simultaneous appearance of several neutrino-produced showers in the field of view of very large 
detectors, such as JEM-EUSO, and 10 TeV gamma radiation from the Virgo cluster. The flux of 
UHE protons from cusps may account for a large fraction of the observed events at the highest 
energies. 
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I. INTRODUCTION 



A. Neutrino astronomy 

Ultra-high-energy (UHE) neutrino astronomy at energies above 10 17 eV is based on new, 
very efficient methods of neutrino detection and on exciting theories for neutrino produc- 
tion. The most interesting range of this astronomy covers tremendously high energies above 
10 19 — 10 20 eV. In fact, this energy scale gives only the low-energy threshold, where the new 
observational methods, such as space-based observations of fluorescent light and radio and 
acoustic methods, start to operate. These methods allow observation of very large areas 
and so detection of tiny fluxes of neutrinos. For example the exposure of the space detector 
JEM-EUSO [1] is planned to reach ~ 10 6 km 2 yr sr. The upper limits obtained by radio 
observations are presented in Fig. [TJ 

The basic idea of detection by EUSO is similar to the fluorescence technique for obser- 
vations of extensive air showers (EAS) from the surface of the Earth. The UHE neutrino 
entering the Earth's atmosphere produces an EAS. A known fraction of its energy, which 
reaches 90%, is radiated in the form of isotropic fluorescent light, which can be detected 
by an optical telescope in space. There is little absorption of up-going photons, so the 
fraction of flux detected is known, and thus EUSOprovides a calorimetric measurement of 
the primary energy. In the JEM-EUSO project [l| a telescope with diameter 2.5 m will 
observe area ~ 10 5 km 2 and will have a threshold for EAS detection E t ^ ~ 1 x 10 19 eV. The 
observations are planned to start in 2012-2013. 

UHE neutrinos may also be very efficiently detected by observations of radio emission by 
neutrino-induced showers in ice or lunar regolith. This method was originally suggested by 
G. Askaryan in the 1960s 0]. Propagating in matter the shower acquires excess negative 
electric charge due to scattering of the matter electrons. The coherent Cerenkov radiation 
of these electrons produces a radio pulse. Recently this method has been confirmed by 
laboratory measurements 0]. Experiments have searched for such radiation from neutrino- 
induced showers in the Greenland and Antarctic ice and in the lunar regolith. In all cases 
the radio emission can be observed only for neutrinos of extremely high energies. Upper 
limits on the flux of these neutrinos have been obtained in the GLUE experiment [4j by 
radiation from the moon, in the FORTE experiment [Bj by radiation from the Greenland 
ice, and in the ANITA j6| and RICE [3] experiments from the Antarctic ice. 

Probably the first proposal for detection of UHE neutrinos with energies higher than 
10 17 eV was made in [8|. It was proposed there to use the horizontal Extensive Air Showers 
(EAS) for neutrino detection. Later this idea was transformed into the Earth-skimming 
effect [9[ for r neutrinos. Recently the Auger detector [Io| put an upper limit on UHE 
neutrino flux using the Earth-skimming effect (see Fig. [1]). 



B. UHE neutrino sources 



What might these new large-area UHE neutrino observatories detect? On the one hand, 
there are without doubt cosmogenic neutrinos, produced by UHECR particles interacting 
with the CMB photons. On the other hand, there may be neutrinos produced in decays or 
annihilation of superheavy particles; this is referred to as the top-down scenario. 



Cosmogenic neutrinos were first discussed in 11]. soon after the prediction of the GZK 



cutoff [12J. There, it was shown that UHE neutrino fluxes much higher than the observed 
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UHECR flux can be produced by protons interacting with CMB photons at large redshifts. 
The predicted flux depends on the cosmological evolution of the sources of UHE protons and 
on the assumed acceleration mechanisms. Recent calculations of cosmogenic neutrino fluxes 
(see e.g. [i3]-(2Q[) are normalized to the observed UHECR flux, with different assumptions 
about the sources. 

The energies of cosmogenic neutrinos are limited by the maximum energy of acceleration, 
E™ x - To provide neutrinos with energies above 1 x 10 20 eV, the energies of accelerated 
protons must exceed 2 x 10 21 eV. For non-relativistic shocks, the maximum energy of ac- 
celeration _E™ ax can optimistically reach 1 x 10 21 eV. For relativistic shocks this energy can 
be somewhat higher. Production of cosmogenic neutrinos with still higher energies depends 
on less-developed ideas, such as acceleration in strong electromagnetic waves, exotic plasma 
mechanisms of acceleration and unipolar induction. 

The top-down scenarios, on the other hand, naturally provide neutrinos with energies 
higher and much higher than 1 x 10 20 eV j2l|. The mechanism common to many models 
assumes the existence of superheavy particles with very large masses up to the GUT scale 



~ 10 16 GeV. Such particles can be produced by Topological Defects (TD) (see |22j for a 
general review). They then rapidly decay and produce a parton cascade, which is terminated 
by production of pions and other hadrons. Neutrinos are produced in hadron decays. 

The production of unstable superheavy particles — the constituent fields of TD — is a 
very common feature of the TD [23J. However, the dynamics of TD is highly nonlinear and 
complicated, the distance between TDs is model-dependent, and the calculation of UHE 



particle fluxes requires special consideration for different types of TD |24| . 

Cosmic strings can release particles in the process of self-interaction, and in the final 
evaporation of tiny loops, but only a few particles are produced by each such interaction. 
Of more interest are cosmic string cusps, where the string doubles back on itself and moves 
with a huge Lorentz factor (25[. Particles emitted by cusps have energies much higher than 



their rest masses, because of the boost. However, the flux from such events is too low to be 



observed 126, 27 



Monopole-antimonopole pairs connected by strings |28l . 1291 . |30| can release superheavy 
particles when the monopole and antimonopole finally annihilate. However, such defects, 
similar to superheavy dark matter (see below), would be accumulated inside galaxies, and 
in particular in the Milky Way. The resulting UHECR flux would be dominated by photons, 
which can reach us easily from short distances. Such photons are not observed [3l[ at 
the level that would be necessary if top-down production were to account for the observed 
UHECR. 

If each monopole is attached to two strings, we have necklaces. Necklaces are an attrac- 
tive source for UHE neutrinos j32j, |33J], but simple models of necklaces may lead to rapid 
annihilation of the monopoles [34|. In other models, however, the monopoles may survive 
for much longer, providing a detectable flux of UHE neutrinos. 1 

In a wide class of particle physics models, cosmic strings can be superconducting, in 
which case they respond to external electromagnetic fields as thin superconducting wires 



1 The main point of Ref. 34[ is that the relativistic motion of strings causes monopoles to develop large 
velocities along the string. As a result monopoles frequently run into one another and annihilate. A 
possible way to avoid this is to consider light strings, which remain overdamped till very late times and 
therefore move slowly. Another possibility is that the strings have zero modes, which act as a one- 
dimensional gas on the strings and slow the monopoles down. These models need further investigation. 
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35l | . String superconductivity arises when a condensate of charged particles (which can be 
either bosons or fermions) is bound to the string. These particles have zero mass in the bound 
state, whereas away from the string they have some mass mx- Loops of superconducting 
string develop electric currents as they oscillate in cosmic magnetic fields. Near a cusp, a 
section of string acquires a large Lorentz boost j c , and simultaneously the string current is 
increased by a factor j c . If the current grows to a critical value J max charge carriers rapidly 
scatter off each other and are ejected from the string. The decay products of these particles 
can then be observed as cosmic rays. This model will be the subject of the present paper. 
Apart from TDs, superheavy particles can naturally be produced by thermal processes 
37| and by time-varying gravitational fields (38l . 39| shortly after the end of inflation. 



These particles can survive until present and produce neutrinos in their decays. Protected by 
symmetry (e.g. discrete gauge symmetry, in particular R-parity in supersymmetric theories), 
these particles can have very long lifetimes exceeding the age of the universe. The resulting 
neutrino flux may exceed the observed flux of UHECR. However, like any other form of 
CDM, superheavy particles accumulate in the Milky Way halo and produce a large flux of 
UHE photons. The non-observation of these photons puts an upper limit on the neutrino 
flux from intergalactic space. 



C. The cascade bound 

The neutrino fluxes are limited from above. The most general upper bound for UHE 
neutrinos, valid for both cosmogenic neutrinos and neutrinos from top-down models, is given 



by the cascade upper limit, first considered in [8|, |40|. The production of neutrinos in these 
scenarios is accompanied by production of high energy photons and electrons. Colliding 
with low-energy target photons, a primary photon or electron produces an electromagnetic 
cascade due to the reactions 7 + 7 ta rget — > e + + e~, e + 7tar g et — > e' + 7', etc. The cascade 
spectrum is very close to the EGRET observations in the range 3 MeV - 100 GeV [4JJ. The 
observed energy density in this range is Regret ~ (2 — 3) x 10 -6 eV/cm 3 . It provides the 
upper limit for the cascade energy density. The upper limit on UHE neutrino flux J u (> E) 
(sum of all flavors) is given by the following chain of inequalities 

An r°° An C°° An 

co cas > — / E'J u {E')dE' > —E / J v {E')dE' = —EJ V {> E), (1) 
c J E c J e c 

which in terms of the neutrino spectrum J V {E) gives 

E 2 J U (E) < —u cas , with u cas < Regret (2) 
An 

Eq. (j2J) gives a rigorous upper limit on the neutrino flux. It is valid for neutrinos produced 
by HE protons, by topological defects, by annihilation and decays of superheavy particles, 
i. e., in all cases when neutrinos are produced through decay of pions and kaons. It holds 
for an arbitrary neutrino spectrum decreasing with energy. If one assumes some specific 
shape of neutrino spectrum, the cascade limit becomes stronger. For a generation spectrum 
proportional to E~ 2 , which is used for analysis of observational data, one obtains a stronger 
upper limit. Given for one neutrino flavour it reads 42| 



E2jt{E) - SArrHEjjE^y (3) 
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FIG. 1: The experimental upper limits on UHE neutrino fluxes in comparison with the electro- 
magnetic cascade upper limit in assumption of E~ 2 generation spectrum (labeled "E 1-2 cascade") 
and with predictions for cosmogenic neutrinos. Neutrino fluxes are given for one neutrino flavour 

Vi + \>i. 



where -Emax and E^n give the range of neutrino energies to which the E~ 2 spectrum extends, 
and i = + or i = v e + v e , or % = u T + u T 

This upper limit is shown in Fig. [IJ One can see that the observations almost reach the 
cascade upper limit and thus almost enter the region of allowed fluxes. 

The most interesting energy range in Fig. [T] corresponds to E v > 10 21 eV, where accel- 
eration cannot provide protons with sufficient energy for production of these neutrinos. At 
present the region of E u > 10 21 eV, and especially E u ^> 10 21 eV is considered as a signature 
of top-down models, which provide these energies quite naturally. 



D. Model assumptions 



In this paper we consider superconducting string loops as a source of UHE neutrinos. 
We consider a simple model in which a magnetic field of magnitude B, occupying a fraction 
of space fs, is generated at some epoch z max ~ 2-3. The strings are characterized by two 
parameters: the fundamental symmetry breaking scale r\ and the critical current J max . We 
take the mass per unit length of string to be fi = rf. 

The predicted flux of UHE neutrinos depends on the typical length of loops produced 
by the strin g ne twork. This issue has been a subject of much recent debate, with different 
simulations [43l . 0, 45, and analytic studies [47, Us) 



yielding different answers. Here 
we shall adopt the picture suggested by the largest and, in our view, the most accurate 



simulations of string evolution performed to date [45|, 4(|. According to this picture, the 



characteristic length of loops formed at cosmic time t is given by the scaling relation 

I ~ at, 



(4) 



with a ~ 0.1. 
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For simplicity and transparency of the formulae obtained in this paper we use several 
simplifications. We assume cosmology without A term with Q c dm + &>b = 1, the Hubble 
constant H = 72 km/s Mpc, age of the universe to = (2/3)H ' 1 , t eq ~ 1 x 10 12 s, and 
(1 + z) 3 / 2 = t /t for the connection of age t and redshift z in the matter era. 

We also assume the fragmentation function for the decay of superheavy X particle into 
hadrons is 

dN/dE oc E~ 2 , (5) 

while Monte Carlo simulation and the DGLAP method give closer to E~ im |49| . 

These simplifications give us a great advantage in understanding the dependence of calcu- 
lated physical quantities on the basic parameters of our model, in particular on fundamental 
string parameter 77. Our aim in this paper is to obtain the order of magnitude of the flux of 
UHE neutrinos and to indicate the signatures of the model. We believe our simplified model 
assumptions are justified, given the uncertainties of string evolution and of the evolution of 
cosmic magnetic fields. 



II. PARTICLE EMISSION FROM SUPERCONDUCTING STRINGS 
A. Particle bursts from cusps 



As first shown by Witten [35(, cosmic strings are superconducting in many elementary- 
particle models. As they oscillate in cosmic magnetic fields, such strings develop electric 
currents. Assuming that the string loop size is smaller than the coherence length of the field 



I < Ib ~ IMpc, the electric current can be estimated as [22j, |35| 

J-O.le 2 ^/. (6) 
Particles are ejected from highly accelerated parts of superconducting strings, called 



cusps, where large electric currents can be induced [50|, |51|- The current near a cusp region 
is boosted as 

J cusp ~ 7 c J, (7) 

where J is the current away from the cusp region and 7 C is the Lorentz factor of the cor- 
responding string segment. Particles are ejected from portions of the string that develop 
Lorentz factors 

7c ~ Jmax/ J, (8) 

where the current reaches the critical value J max- This maximum current is model- 
dependent, but is bounded by J ma x ^ er], where r\ is the symmetry breaking scale of the 
string and e ~ 0.1 is the elementary electric charge, renormalized to take into account 



self-inductance 22 



One may parametrize J maa; by introducing the parameter i c < 1: 

Jmax ^c&Vi (9) 

If the charge carrier is a superheavy particle X with mass nix, the case which will be 
considered here, one may use e r x for the energy of X-particle in the rest system of the cusp 
and ex in the laboratory system. Then e r x = ^mx = i c V and 

ex ~ idcV, (10) 
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respectively, where 7 is the average Lorentz factor of X-particle in the rest system of the 
cusp. In Eq. (|T0|) we took into account that that the energy of X-particle in the laboratory 
system is boosted by the Lorentz factor of the cusp j c . 

The number of X particles per unit invariant length of the string is ~ J/e, and the 
segment that develops Lorentz factor 7 C includes a fraction 1/7 C of the total invariant length 
I of the loop. Hence, the number of X particles ejected in one cusp event (burst) is 

N b x ~ (J/e)(// 7c ) ~ J 2 l/eJ max . (11) 

The oscillation period of the loop is 1/2, so assuming one cusp per oscillation, the average 
number of X particles emitted per unit time is 



N x ~ 2J 2 /eJ max , (12) 



and the luminosity of the loop is 



->tot 



N x e x - (13) 



The X particles are short-lived. They decay producing the parton cascade which is 
developed due to parton splitting in the perturbative regime, until at the confinement radius 
the partons are converted into hadrons, mostly pions and kaons, which then decay producing 
gamma rays, neutrinos, and electrons. These particles together with less numerous nucleons 
give the observational signatures of superconducting cusps. 

The neutrino spectrum at present epoch z = 0, produced by the decay of one X-particle 
with energy e x ~ iclcV at epoch z can be calculated using the fragmentation function (J5]) 
for an X-particle at rest: 

where £™* and E™** are the maximum and minimum neutrino energies in the rest system 
of X-particle. 

Particle emission from a cusp occurs within a narrow cone of opening angle 

Oc ~ % X ~ J/Jmax (15) 



The duration of a cusp event is [51 



Uurst ~ he 3 ( 16 ) 



B. Superconducting loops in the universe 

In any horizon-size volume of the universe at arbitrary time there are a few long strings 
crossing the volume and a large number of small closed loops. As loops oscillate under the 
force of string tension, they lose energy by emitting gravitational waves at the rate 

E g ~ TGfj 2 , (17) 

where /1 ~ rj 2 is the string mass per unit length, G = l/m 2 Pl is the gravitational constant 
and T ~ 50 is a numerical coefficient. 
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The number density of loops with lengths in the interval from I to I + dl at time t can 
be expressed as n(l,t)dl. Of greatest interest to us are the loops that formed during the 
radiation era t eq and still survive at t > t eq . The density of such loops at time t is given by 

n(/,t)rf/~^ 2 r 2 /- 5/2 rf/, (18) 
in the range from the minimum length l min to the maximum length I ~ at eq , where 

U ~ TGfjt ~ 3 x 10 n r/ 2 (l + ^ 3/2 cm (19) 

and r] 10 = ?7/10 10 GeV. Here and below we assume that the loop length parameter in (TJJ 
is a ~ 0.1, as suggested by simulations j45|, 46[. Loops of the minimum length are of most 
importance in our calculations because they are the most numerous. 

For a loop of length I at redshift z, the Lorentz factor at the cusp -y c can be expressed as 

mm /r>n\ 

lc = J = ^Bl =lc{lmm) — (20) 



where y c (l min ) = y (l + z) 3/2 and 



10i c ri 6 x 10 n z c , . 

eBt TG/j, B_ 6 t] 10 



where B_ 6 is the magnetic field in microgauss. 



C. Limits on rj 

The string motion is overdamped at early cosmic times, as a result of friction due to 
particle scattering on moving strings. The friction-dominated epoch ends at 

U ~ (Gfi)-% (22) 

where t p is the Planck time. In the above analysis we have assumed that loops of interest 
to us are formed at t > £*. The corresponding condition, 

TGfito/a > U, (23) 

yields 

r] > 10 9 GeV. (24) 
For strings with rj < 10 9 GeV we should use 

lmin ~ at* (25) 

instead of (fl9|) for the length of the smallest surviving loops. 

We should also verify that energy losses due to particle emission and to electromagnetic 
radiation in recent epochs (after magnetic fields have been generated) are sufficiently small, 
so the lifetimes of the loops (which we estimated assuming that gravitational radiation is 
the dominant energy loss mechanism) are not significantly modified. 

The average rate of energy loss due to particle emission is 

E part ~ 2f B N x e x ~ f B JJma X /e 2 (26) 
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where we have used Eqs. fTTTT) and (jTOl) . The electromagnetic radiation power is smaller by 
a factor e 2 ~ 10 -2 . 

The factor f B in Eq. (|26|) is the filling factor - the fraction of space filled with the magnetic 
field. It gives the fraction of time that cosmic string loops spend in magnetized regions. We 
assume that loop velocities are sufficiently high that they do not get captured in magnetized 
cosmic structures (such as galaxy clusters or LSS filaments). To justify this assumption, we 
note that particle emission can start only after the cosmic magnetic fields are generated, 
that is, at z ~ 3 or so. Before that, gravitational radiation is the dominant energy loss 
mechanism, and the loops are accelerated to high speeds by the gravitational rocket effect 
I52I . 53]. The smallest loops of length ffl9|) have velocities v ~ 0.1, certainly large enough to 



avoid capture. 

The particle emission energy rate (1261) should be compared to the gravitational radiation 
rate 

E g ~ FG/j 2 . (27) 

The ratio of the two rates is zero at z > z max , where z max ~ 2-3 is the red-shift of magnetic 
field production. At z < z max it is given by 

E part /E g ~ 25iL 3 £U^r 1 (t-) (! + z )' 3/2 - (28) 

\ <"min J 

where /_ 3 = /s/10 -3 and l min is given by ([19]). 

If particle emission is the dominant energy loss mechanism, then the lifetime of a loop is 

$ IO77 t r]w , QO s 

E pa rt ez c f B B f- 3 B- 6 i c 

Note that r is independent of I. This means that all loops surviving from the radiation era 
decay at about the same time. 

For the time being, we shall assume that particle radiation is subdominant. We shall 
discuss the opposite regime in Section II. G. 



D. Rate of cusp events 

The rate of observable cusp bursts (i.e., the bursts whose spot hits the Earth) is given by 

. dQ . --dV(z) 
dN b = f B — z)dl—y- 30 

47T 1 + Z 

where, as before, fg is the fraction of space with magnetic field B, dQ = l-nQdB is the solid 
angle element, with 9 limited by the angle of cusp emission 9 C ~ 1/7 C ; z) = n(l, z) / (1/2) 
is the frequency of the bursts with n(l, z) given by Eq. ffl8|) . and dV(z) is a proper volume 
of space limited by redshifts z and z + dz, 

dV(z) = 547rfj|[(l + z) 1 ' 2 - 1] 2 (1 + z)- u ' 2 dz. (31) 

Integrating Eq. (|30|) over 9, I and z, we obtain 

N b = 54 7 r(t e(? t ) 1 / 2 (rG yU )- 1 / 2 (e/10z c r ? ) 2 dz [{1 ^^'^ f B (z)B 2 (z), (32) 
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where z max is the redshift at which the magnetic fields are generated. Since the earth 
is opaque to neutrinos with the energies we are considering, only half of these bursts can 
actually be detected by any given detector at the surface of the earth or using the atmosphere. 

The value of the integral in (j3"2~j) depends on one's assumptions about the evolution of 
the magnetic field B and of the volume fraction fg. This evolution is not well understood. 
If we take these values out of the integral in Eq. (1321) as the average and characterize them 
by the effective values of parameters B_ 6 and /_3 in the range < z < z max , then Eq. 
reduces to 

N b = 2.6 x 10 2 ".- 6 3~ 3 n ^ yr-\ (33) 



2 B 2 _ 6 f_ 3 I 

where the integral 



'rfio 0.066 



W ^ [( \t+I)V ]2 = |[l-(l+-)- 3/4 ]-|[l-(l+^)- 5/4 ] + ^l-(l+-)~ 7/4 ] ? (34) 

is equal to 0.015, 0.042 and 0.066 for z = z max — 1, 2 and 3, respectively. 

The integrand in Eq. f[3"2"l) includes the product f B (z)B 2 (z). In the calculations of other 
physical quantities below, similar integrals will have different combinations of fs{ z ) an d 
B(z). Nevertheless, we shall assume that the average values taken out of the integral are 
characterized by approximately the same values of fs and -B_6- 

All cosmic structures — galaxies, clusters, and filaments of the large-scale structure - 
are magnetized and contribute to the rate of cusp bursts. In the recent epoch, z < 1, the 
dominant contribution is given by clusters of galaxies with i?^ 6 /_3 ~ 1. The magnetic fields 
of galaxies have about the same magnitude, but the corresponding filling factor f B is orders 
of magnitude smaller. We shall assume that this holds in the entire interval < z < z max . 
The sources in our model are then essentially clusters of galaxies. 



E. Diffuse flux of UHE neutrinos 

Including all produced neutrinos, the diffuse differential neutrino flux can be given by the 
formula 

M E) = lJ d N b N> x UE) 1 ^ (ry (35) 

where dN b is the rate of neutrino bursts (13"U|) . is the number of X particles produced 
per burst, given by Eq. (fTTj) . £„(E) is the neutrino spectrum produced by the decay of one 
X-particle, given by ( 1T41) . 

Q jet = n6 2 c = ^, (36) 

r(^) = 3t [l-(l + ^r 1/2 ] (37) 

is the distance between a source at redshift z and the observation point at z = 0, and Qj et r 2 
is the area of the burst spot at the Earth from a source at redshift z. 
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Using expressions ([TBI) and (I3TT) . and assuming that the product fB(z)B(z) does not 
change much in the interval < z < z max , we obtain 2 



F 2 t m _ 0.3i c m p i(t eq /t )(eBt Q )f B 7/4 



nun J 



Numerically, this gives for the neutrino flux summed over neutrino flavors 

E 2 J U (E) = 6.3 x l(T 8 i CJ B_ 6 ,/l3 Gel/ cm" 2 s' 1 sr -1 , (39) 

where we have set z max = 3 and estimated the logarithmic factor as ~ 30. 

For i c ~ 1, the flux (1391) is close to the cascade upper limit shown in Figure [TJ Notice 
that the diffuse neutrino flux (1381) does not depend on 77. The neutrino flux must correlate 
with clusters of galaxies. 

To detect this flux, we need to monitor a target with some large mass M. The effective 
cross-section of the detector is then 

£ = a vN M/m N (40) 

where a V N ~ 3 x 10~ 32 cm 2 is the neutrino-nucleon cross section at E > 10 10 GeV and 
the mass of a nucleon. Because of the opacity of the earth, the detector will see solid angle 
about 2tt sr. The detection rate of particles with energy above E is 

2itEJ v {E)Y> = 21 ( ^GeV ^ -^ ^ ^-i (41) 

In the case of JEM-EUSO in tilt mode, M ~ 5 x 10 18 g, and thus we expect about 100i c 
detections per year, so events can be expected for i c > 0.01. 



F. Neutrino fluence and the number of neutrinos detected from a burst 

The fluence of neutrinos incident on the detector from a burst at redshift z can be 
calculated as 

*(> E) = f"> f (42) 

Consider a neutrino burst from a loop of length / at redshift z. Using N b x from ( flTi) . l min 
from ( Ti~9i) and ^(> E) from ( TT4T) . we obtain for a loop of any length /, 

$(> E) = m ^ (43) 

1 > lMeBtlEHE^/E^)[{l + z)^-ir 1 ' 

which numerically results in 



2 We note that numerical simulations of the magnetic field evolution performed by Ryu et al. [54[ do indicate 
that the space average of the magnetic field (B(z)) = fB(z)B(z) remains roughly constant at ~ 10 -9 G 
for < z < 3 and decreases at larger values of z. The effective values B_§ and /_3 could be different 
from those in Eq. ([3"3")1 for the rate of bursts, but we neglect the possible difference, 
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The number of neutrinos detected in a burst is 

Nf* ~ $(> E)E (45) 

With M ~ 5 x 10 18 5f as above, 



N: 



^q.O.U^ ' ,„ ,46, 



£ B_ 6 [(l + z )i/2_i] 



Therefore, for a certain range of i c ?7io values and source redshifts z, multiple neutrinos can 
be detected as parallel tracks from a single burst. For example, for i c r] w ~ 3, and z ~ 1, 
A^ et ~ 17. 

For neutrino energies of interest, E v > 1 x 10 20 eV, the neutrino Lorentz factor is so large 
that there is practically no arrival delay for neutrinos with smaller energies. All neutrinos 
from a burst arrive simultaneously and produce atmospheric showers with parallel axes, 
separated by large distances. 

For other sets of parameters N^ et < 1 , i.e. only one neutrino from a burst (or no 
neutrino) is detectable. As rj increases, the rate of bursts (133]) diminishes while the number 
of neutrinos per burst increases, so that the total neutrino flux remains unchanged. The 
rate of detected neutrino bursts with the number of detected neutrinos N^ et > ( for each 
burst, is given by Eqs (1331 and (1341) . with z max determined by N^ et (> E, z max ) = ( . Using 
Eq. (g6]) we obtain for x max = (1 + z max ): 



1 + 



/ O.ingr^o 10 10 GeV \ 

V C ~b7s e ) 



1/2" 1 2 



(47) 



if ( 1471) is less than 4, and x max = 4 if (j47|) is larger than 4. Introducing in Eq. (I3"3l coefficient 
1/2 which approximately takes into account the absorption of UHE neutrinos crossing the 
Earth we obtain for the rate of detected bursts with N^ et > £ 

Nt l (> = 2-1 x ltfL^I( Zmax ) yr-\ (48) 



where I{z max ) is given by Eq. (|34|) with z max from Eq. (|47|) . 

In Fig. [5], we have shaded the region of the parameter space (rj, i c ) corresponding to a 
detectable flux of neutrinos. Curved lines in the figure mark the regions where we expect 
a burst with a given multiplicity of neutrinos, ( = 2, 3 or 10, detected simultaneously by 
a detector with parameters of JEM-EUSO tilted. To the left of the 2-neutrino-burst line, 
only a diffuse flux of single neutrinos can be observed. This flux depends only on i c , and 
the vertical left boundary of the shaded region marks the value of i c at which it drops below 
one particle detected per year. 

Note that the regions shown for multiple events are those where we expect at least one 
burst per year whose average multiplicity is the given ( or more. But it is possible even if 
the parameters are to the left of the ( = 2 line that we would happen to observe multiple 
neutrinos from a single burst, which would give a clear signature of neutrino-jet emission 
from cusps. 

Another quantity of interest is the rate of detected neutrinos f v (> C) i n the events with 
neutrino multiplicity greater than (. It is given by 
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0.01 0.1 1 



'c 

FIG. 2: The region of parameter space where neutrinos can be seen by a detector with the param- 
eters of JEM-EUSO. The curved lines show the left edges of the regions in which bursts containing 
at least 2, 3, and 10 neutrinos can be expected at least once per year. Below the dotted line, 
particle radiation is the dominant channel of energy loss from loops. 

The important feature of the calculations is the independence of iV^ e *(> E, z, I) from I. This 
allows us to integrate over I in Eq. fH9|) to obtain 

/„(> C) = 2.1 x 10 3 ^^ / dz(l + z)~ [(1 + z) 1 ' 2 - 1] N d v et {> E, z), (50) 

l c r llQ Jo 

where z max (() is given by Eq. fH7|) . Using Eq. (H6T) for N^ et (> E,z) results in 

/„(> = 1-3 x 10\f_ 3 B_ 6 [l - x^l\i c , Vl0 )} yr" 1 . (51) 
The asymptotic expression at 0.11z 3 r/f /i?_ 6 C <C 1 gives 

U> = '-^^^B^ yr-. (52) 

G. Neutrino fluxes in the particle-emission dominated regime 

So far we have assumed that gravitational radiation is the dominant energy loss mecha- 
nism of strings. In the opposite regime, where the particle emission energy losses dominate, 
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the loop's lifetime T part is independent of its length and is given by Eq. (1291) , 



10r? 

eicjBB 

We shall analyze this regime in the present section. 

As before, we shall adopt the idealized model where the magnetic field B is turned on at 
time t = ts, corresponding to redshift z max , 

t B ~t Q (l + z max )- 3/2 . (54) 

The loops decay at the time tdec ~ + T part . The rate of observable bursts iV& is given by 
Eq. (133]) with I from Eq. AMI) , where the integration is taken between z dec and z max and z dec 
is the redshift corresponding to the time t dec . 

If r part > the redshift z dec is significantly different from z max , with = z max - z dec > 
1, and the value of I is not much different from that evaluated in Sec. II. D. This is an 
intermediate regime, in which the results we obtained in Sections II. D and II. E for the rate 
of bursts and for the diffuse flux can still be used as order of magnitude estimates. 

For T part <C ts, the loops lose all their energy to particle emission in less than a Hubble 
time. The condition T par t ~ ts can also be expressed as E par t/ E g (z max ) ~ 1, or, using 
Eq. ([28]) with z max ~ 3, 

7] ~ 3 x 10 10 2 c /_ 3 5- 6 GeV. (55) 

It marks the boundary of the strong particle-emission domination regime and is shown by 
the inclined dotted line in Fig. 2. Below this line, the results of the preceding sections do 
not apply even by order of magnitude, but as we shall see, detectable neutrino fluxes can 
still be produced. 

The redshift interval Az = z max — z dec for T par t <C £b can be estimated as 

As«?^(l + w)<l, (56) 
and the integral / in Eq. (|34"|) is given by 

j — A . z max ) ij ( . 

(1 + Z ) n /4 • ^ 

With z max ~ 3, we have ts ~ to/8, and 

lESH „ 0,4 - V '° , . (58) 

tB J-3-D-6^c 

The rate of bursts that are actually detected, N^ 1 , can be expressed as a product of N b 
and the probability pf, et that at least one neutrino from the burst will be detected. This 
probability is simply related to the average number of detected neutrinos per burst N^ et , 
given by Eq. f|46|) . 

= l-exp(-iV^). (59) 

For N^ et < 1, we have 

p det ^ N det ^ 



14 



and 

Nft ~ NtNf* ~ ]_ 20r]l \ 7/4 yr- 1 ~ lOryxo yr' 1 , (61) 

I, ' T Zmax ) 

where in the last step we have used z max ~ 3. Requiring that N^ et > 1 yr^ 1 , we obtain the 
condition 

T] > 10 9 GeV. (62) 

Note that at the boundary of detectability, where r] ~ 10 9 GeV, we always have N^ et <C 1, 
and thus the approximation (I60I) is justified. This boundary is the lower horizontal line 
bounding the observable parameter range in Fig. 2. Note also that Eq. (16"2l coincides with 
with the condition (1241) for the burst-producing loops to be unaffected by friction. 

It is interesting to note that the detection rate (EH) in the particle-emission dominated 
regime is independent of i c and depends only on the symmetry breaking scale rj. This is in 
contrast with Eq. (T4TT) for the case of gravitational radiation dominance, where the rate is 
proportional to i c and independent of rj. 



H. Cascade upper limit on neutrino flux in the superconducting string model 

In Subsection I.C of the Introduction, we gave a very general upper limit for UHE neutrino 
flux. The presence of such a limit does not contradict the existence of stronger upper limits 
in some particular models with additional assumptions. 

In this section, we calculate the energy density of the cascade radiation in our model and 
compare it with uj cas = 2 x 10~ 6 eV cm~ 3 allowed by EGRET measurements. 

The cascade energy density can be calculated as 

/dz r l max(z) 
—— — rj / dlf B n(l,t)L em (l,t) (63) 
\ L + Z > Jl mi „(z) 

where L em (l, t) = ^L tot (l, t) is the loop luminosity in the form of UHE electrons and photons 
produced by pion decays. The standard calculation results in 

Ucas = 7(TGfiy/ 2 4 = l cf-sB_ 6 eVcm (64) 

The energy density AMI) does not depend on rj and since u cas < uj egret , it respects the 
general upper limit ([3]). For i c ~ 1, the predicted neutrino flux (13"9|) is close to the upper 
limit shown in figure [TJ 



III. GAMMA-RAY JETS AND SINGLE GAMMA-RAYS FROM THE CUSPS 
A. Bursts from loops in the Milky Way 

In each galaxy, including the Milky Way, there are approximately Ni loops with I > l m i n , 

N t ~ n{> l min )V g ~ 2.5 x lOVo 3 Vl0 3 fcpc 3 , ( 65 ) 

where V g is the volume of the magnetized part of the galaxy. A narrow jet of particles 
emanating from a cusp on such a loop can in principle hit the Earth. The probability of 
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such a catastrophic event is very small because of the smallness of solid angle Vtj et of jet 
emission. The number of jets hitting an area S on the Earth per unit time does not depend 
on S if S <C Qj e tr 2 , where r is the distance from the source. This rate is given by 

N b = PV g jdl^fl, (66) 

where again we assume one cusp per oscillation, and P = Qj e t/47t = l/(47 c ) 2 is the proba- 
bility to hit the detector. After the standard calculations, we obtain for V g ~ 1 x 10 3 kpc 3 , 

N b = lx 10- 13 B 2 _ 6V ^z- 2 yr' 1 - (67) 

Thus, for particles propagating rectilinearly, the jets from cusps in our galaxy are unob- 
servable. 

The most important components of the galactic jets are photons and neutrinos. A photon 
jet at the highest energies undergoes widening of the jet angle due to photon absorption in 
the galactic magnetic field [55|. Absorption of photons, 7 + B —> B + e + + e~, is followed by 
energy loss by electrons and positrons in the magnetic field, with the emission of synchrotron 
photons in directions different from that of the primary photon. This results in the widening 
of the solid angle Qj et (55[. 

The widening of photon jets in the Milky Way is negligible. This can be illustrated by a 
numerical example. The highest energy of a photon in a jet is E™ ax ~ j c t] ~ 10 31 i c /B_ G eV. 
Photons with E 1 ^ 10 25 eV are absorbed in galactic magnetic fields. The produced electrons 
and positrons with E e ~ 10 25 eV have lifetime r ~ 10 3 s for synchrotron energy losses 
and attenuation length l att ~ 3 x 10 13 cm. Since the Larmor radius of such electrons is 
r L ~ 3 x 10 28 cm, the deflection angle 6 ~ l a tt/ r L ~ 10~ 15 is of no consequence. 



B. Cascade gamma- radiation from Virgo cluster 

As was discussed above, the photon jets from the galactic cusps are not widening and 
thus are invisible. For cusps at large distances, the widening of the photon jet efficiently 
occurs in the cascading on CMB photons, 7 + Jcmb — > e + + e~, e + Jcmb — > e' + 7' etc., 
and the source can be seen in gamma-radiation. As in the case of diffuse cascade radiation 
(see section I.C), all primary photons with energy higher than the absorption energy e a are 
absorbed on CMB radiation and converted into low-energy cascade photons. Thus, cusps 
can be seen in 100 GeV — 100 TeV gamma radiation, similar to the sources of UHE protons 
which can be seen in TeV gamma-radiation [56j. 

The nearest source from which this radiation can be expected is the Virgo cluster. It 
is located at distance r = 18 Mpc, and the number of loops within the core of radius 
R c ~ 3 Mpc, where a magnetic field B ~ 10~ 6 G can be reliably assumed, can reach 
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niR c ~ 7 x 10 rj 10 , , with the luminosity in e + e component for each loop L] oop ~ 2 
10 29 i c r/f B~ 6 ergs^ 1 . Half of this energy goes into cascade radiation, L cas ~ 0.5L] oop . 
The s pec trum of the cascade photons at distance r ~ 20 Mpc has two characteristic 



energies |40jj: the absorption energy e a ~ 100 TeV and the energy e x . The latter is the 
energy of a photon produced by an electron born in 7 + Jcmb — ► e+ + e ~ scattering by a 
photon with = e a . The energy of this electron is E e ~ 0.5e a ~ 50 TeV and the second 
characteristic energy is e x ~ 20 TeV for r ~ 20 Mpc. 

The spectrum of cascade photons at observation is calculated in [401 as 
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The spectrum constant K in (1681) can be expressed in terms of the cascade luminosity 
L cas and the distance to the source r as 



K = Lcas - (69) 

where £l e ff is the effective solid angle produced by scattering of cascade electron in extra- 
galactic magnetic field. In case of full isotropization Q e ff ~ 47r. Cascade luminosity can 
be estimated as 1/4 of the total luminosity of cusps in a cluster, L cas ~ \Li oop Ni oop . Using 
Lioop — 4.4 x 10 29 ?7^ i5_ 6 erg s^ 1 and Ni oop ~ 2.5 x lO 11 ?]^ 3 , valid for a cluster core with 
R c ~ 3 Mpc, one obtains for the flux 

J 7 (> 6x ) = / rfE 7 J 7 (E 7 ) ~ 1 x l(T 13 z CJ B_ 6 (i? c /3 Mpcf cm- 2 s' 1 (70) 

J £i 

which is marginally detectable by present telescopes. Note that L cas and the flux J 7 do not 
depend on 77. We consider the estimate ATP]) as a very rough indication of detectability of 
the gamma-ray flux from the Virgo cluster. Much more accurate calculations are needed for 
a reliable prediction of this flux. 



IV. UHE PROTONS FROM SUPERCONDUCTING STRINGS 

The cusps of superconducting strings in clusters of galaxies produce UHE nucleons at 
fragmentation of parton jets with a fraction of nucleons €n = 0.12 [33[ relative to the total 
number of hadrons. The generation rate Q P {T P ) of UHE protons with Lorentz factor T p per 
unit comoving volume and unit time can be expressed through emissivity, 

■pmax 
I P 

C = dT p m N T p Q p (T p ) , (71) 

/ rmin 
P 

where the emissivity Co is the energy released in UHE protons at z = per unit comoving 
volume per unit time, r™ a3: and T p nm ~ 1 are the maximum and minimum Lorentz factors of 
the protons, respectively, and is the nucleon mass. For a power-law generation spectrum 
Q P {T P ) ~ T" 2 , we have 



The emissivity is calculated as 

dln(l)L cusp (l), (73) 

- iTiin 

where l min is given by (fl9j) . while n(l) and L cusp are given by ffT8]) and fflBl respectively. For 
L cusp one readily obtains 
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Lcusp = JJYj^T = ^.2i c eBlr], (74) 
and after a simple calculation we have 

Co = 0Az c e N f B {te PJ^" J- J^- 3 = 1-4 * W*\f„ 3 B_ 6 erg Mpc' 3 yr~\ (75) 

One more parameter relevant for the calculation of Q P (T P ) is T™ ax = E™ ax /mN. It can 
be estimated using E™ ax ~ O.lex, where ex = icIcV is the energy of the boosted X particles 
in the laboratory system, which being estimated for loops of length l m i n , gives 



r r = 1 x (76) 



Notice that r™ aa: does not depend on rj and that it enters Q P (T P ) through lnr™ aa; . 

Now we can calculate the space density of UHE protons using the generation rate Q P (T P ) 
given by (1721) and taking into account propagation through CMB radiation with the help of 
the kinetic equation 

d d 

— n p (T p , t) - — [b(T p , t)n p (T p , t)] = Q p (T p , t), (77) 
at o Tp 

where b(T p ,t) = —dT/dt describes energy losses of UHE protons interacting with CMB 
photons. For T ^ 3 x 10 10 , the proton energy losses become large and one can neglect the 
first term in the lhs of equation (1771) . Then Eq. (1771) becomes stationary and its solution for 
t = t reads 

"* (ry = 6(ry/ rp «' (r ' ),ff ' = mw r,4(r°)taiw (78) 

In terms of the proton energy E = m^Tp and the diffuse flux J P (E) = (c/4-K)n p (E), we 
have, in the standard form of presentation, 

E3jp{E) = £inr^r p 6(r p )- (79) 

A numerical estimate at E = 3 x 10 19 eV gives 

E 3 J P (E) = 1.3 x 10 24 z c /_ 3 5_ 6 eV 2 m~ 2 s' 1 sr' 1 . (80) 

With i c ~ 1, the calculated flux (|80|) coincides well with the measurements at the same 
energy, e.g., with the HiRes flux E 3 J P (E) = 2. Ox 10 24 eV 2 m~ 2 s' 1 sr~ x , so the cusp emission 
may account for the observed events at the highest energies. For i c < 0.1 the UHE proton 
flux is considerably less than the observed one. 

The UHE proton spectrum from superconducting strings has a sharper GZK cutoff than 
the standard spectrum for homogeneously distributed sources. This is due to the absence 
of clusters of galaxies in the vicinity of our galaxy. The nearest cluster, Virgo, is located at 
18 Mpc from the Milky Way; other clusters are located at much larger distances. Nearby 
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sources affect the spectrum at E ^ 1 x 10 20 eV, where the proton spectrum from supercon- 
ducting strings is predicted to be steeper than the standard one. The experimental data at 
present have too low statistics to distinguish the two cases. 

In contrast, homogeneously distributed sources such as necklaces [33 |. give the dominant 
contribution at E ^ (7 — 8) x 10 19 eV in the form of UHE photons, coming from nearby 
sources. In the case of superconducting strings such component is absent. The UHE photon 
component from superconducting strings is not dominant at energy lower than 5 x 10 19 eV, 
because absorption of photons at these energies is stronger than for protons. 



V. CONCLUSIONS 

Superconducting cosmic strings produce high energy particles in the decay of charge 
carriers, X particles, ejected from the string cusps. The large Lorentz factor 7 C of the cusp 
boosts the energies of these particles and collimates them in a narrow beam with opening 
angle 9 ~ l/7c- The basic string parameter is rj, the scale of symmetry breaking, which 
we parametrize as rj = ?7iol0 10 GeV. Another free parameter i c < 1 determines the critical 
electric current in the cusp, J ma x = ic^Vi an d the mean energy of the charge carriers X 
escaping from the string, Ex = i c V- 

The astrophysical parameter which determines the electric current induced in the string 
is the magnitude of the magnetic field B in the relevant cosmic structures. The fraction 
of the universe occupied by magnetic field B determines the flux of high-energy particles 
produced by superconducting strings. The most favorable values of B and /# for the gen- 
eration of a large flux of UHE neutrinos are B ~ 10~ 6 G and /# ~ 10~ 3 . They correspond 
to clusters of galaxies. 

The main uncertainties of our model are related to the uncertainties in our understanding 
of the evolution of cosmic strings and of the origin and evolution of cosmic magnetic fields. 
On the cosmic string side, the key unknown quantity is the parameter a which sets the 
characteristic length of string loops in Eq. (T4J). Here, we used the value of a ~ 0.1, as 
suggested by numerical simulations in Refs. [45l . We have also disregarded the effects of 



loop fragmentation. Toward the end of its life, the loop's configuration may be sufficiently 
modified by radiation back-reaction that the loop will self-intersect and break up into smaller 
pieces. These smaller loops will have more frequent cusps, shorter lifetimes, higher velocities, 
and smaller induced currents. The effect of such loops on the neutrino fluxes is hard to 
assess without a quantitative model of loop fragmentation. This will have to await the next 
generation of string evolution simulations. For now, our results in this paper can be regarded 
as lower bounds on the neutrino fluxes from superconducting strings. 

On the astrophysical side, basically unknown is the cosmological evolution of the magnetic 
field parameters /b(V) and B(z) in the redshift interval < z < z max , where z max ~ 2 - 3 is 
the redshift when the magnetic field was generated. For the space average value (fB(z)B(z)) 



we use the numerical simulation by Ryu et al. [541 ]. according to which this value remains 
roughly constant at < z < 3. Some important quantities, such as the diffuse neutrino 
flux J„(E), the cascade energy density u cas , and the UHE proton emissivity are determined 
by the evolution of the product fs(z)B(z). However, some other quantities, such as the 
rate of neutrino bursts and fluence depend on the evolution of /b(z) and B(z) in other 
combinations. In these cases we consider the parameters /_3 and B_ e as effective values, 
using /_ 3 ~ B_ G ~ 1. 

In addition, we adopted the following simplifying assumptions. The Lorentz factor of the 
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cusp is characterized by a single fixed value 7 C , while in reality there is a distribution of 
Lorentz factors along the cusp. The spectrum of particles in a jet is approximated as E~ 2 , 
while a QCD calculation j49| gives a spectrum which is not a power law, with the best power- 
law fit as E~ im . We use cosmology with A = 0. The spectrum of photons from Virgo cluster 
and the diffuse spectrum of UHE protons are calculated using very rough approximations. 
Given the uncertainties of string and magnetic field evolution, these simplifications are 
rather benign. On the other hand, they have the advantage of yielding analytic formulae, 
which allow us to clearly see the dependence of the results on the parameters involved in 
the problem. In particular, with the assumed particle spectrum ~ E~ 2 , the diffuse flux 
of neutrinos, the cascade upper limit, the flux of TeV photons from Virgo cluster and the 
diffuse flux of UHE protons do not depend on rj. Since the realistic spectrum is very close 
to E~ 2 , this means that the quantities listed above depend on rj very weakly. 
We summarize the results obtained in this work as follows. 

As our calculations show, among different sources, such as galaxies, group of galaxies, 
filaments, etc., the largest diffuse flux is produced by clusters of galaxies with B ~ 10~ 6 G 
in a cluster core and fs ~ 10~ 3 . The calculated diffuse neutrino flux for three neutrino 
flavours is 

E 2 J U (E) ~ 6.3 x KrtjLgS-e GeV cm~ 2 s' 1 sr -1 . (81) 

This flux respects the cascade upper limit, provided by the energy density of electrons, 
positrons and photons, which initiate electromagnetic cascades in collisions with CMB pho- 
tons. The cascade energy density is calculated as 

u cas = 9.5 x lQ-\f_ z B 2 _ % eVcm- 3 . (82) 

and is close to the cascade limit for i c ~ 1. 

At energies E < 10 22 eV, the flux flS]) is detectable by future detectors JEM-EUSO 
and Auger (South + North). The signature of the superconducting string model is the 
correlation of neutrinos with clusters of galaxies. We note, however, that the neutrino flux 
from the nearest cluster, Virgo, is undetectable. 

Another signature of the model is the possibility of multiple events, when several showers 
appear simultaneously in the field of view of the detector, e.g. JEM-EUSO. They are 
produced by neutrinos from the same jet. The time delay in arrival of neutrinos with 
different energies is negligibly small. Such multiple events are expected to appear for a 
certain range of parameters, as indicated in Fig. 2. 

As an illustration, in Tabled we show, for a representative value r\ = 5 x 10 10 GeV, the 
diffuse neutrino flux, in units of the cascade upper limit J™ ax , the rate of bursts, and the 
average shower multiplicity for several values of i c . Note that the bottom row in the table 
is the average multiplicity, that is, the average number of neutrinos detected per burst. For 
example, the low multiplicity at i c = 0.1 indicates that only a small number (about 3) out 
of the 140 bursts per year will actually be detected. For i c = 1/3, the average multiplicity 
is below 1, but Fig. 2 shows that we can expect at least one 2-neutrino burst per year. 

A photon jet from the cusp initially propagates together with the neutrino jet, within the 
same solid angle. However, at large enough distance, photons from the jet can be absorbed 
in collisions with CMB photon (7 + jcmb e + + e~), the produced electrons (positrons) 
emit high-energy photons in inverse- Compton scattering (e + jcmb — » e' + 7'), and thus 
an em cascade develops. Electrons are deflected in magnetic fields, and photon radiation is 
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TABLE I: The diffuse flux J U (E) in units of the cascade upper limit J™ ax , the rate of neutrino 
bursts, and the shower multiplicity (the average number of neutrinos detected in one bursts), 
for 7] = 5 x 10 10 GeV and different values of i c . The multiplicity is shown for neutrinos with 
E > 10 10 GeV from a burst at z = 2. 



i c 


1.0 


1/2 


1/3 


0.1 


JjJ? ax 


0.38 


0.19 


0.13 


0.038 


rate of bursts 


1.4 yr -1 


5.4 yr -1 


12.2 yr" 1 


140 yr" 1 


multiplicity 


21 


11 


0.8 


0.022 



isotropized. Due to this effect, 10 — 100 TeV gamma radiation from the nearby cluster of 
galaxies, Virgo, can be marginally detectable. The corresponding photon flux is given by 

J y (> ej~lx 10~ 13 vB_ 6 cm- 2 s- 1 (83) 

where e x ~ 20 TeV. 

In the Milky Way, there may be a large number of loops with radiating cusps, but because 
of the very small jet opening angle, the probability to observe UHE particle jets coming from 
these loops is extremely small. 

The diffuse flux of UHE protons is suppressed by the small fraction of nucleons produced 
at decay of X particles (the factor = 0.12 is obtained in MC and DGLAP calculations 
|49j). and by energy losses of protons interacting with the CMB during propagation. The 
calculated flux at energy E ^ 3 x 10 19 eV is given by the approximate formula 

E3 ME) - £i^p^y (84) 

where b(E) = —dE/dt is the energy loss rate of protons, r™ ax is the maximum Lorentz 
factor of a proton at production, and Cq is the emissivity (energy in the form of protons 
emitted per unit comoving volume per unit time), given by 

C = 1.4 x 10 45 *J_ 3 5_ 6 erg Mpc' 3 yr' 1 (85) 

For i c ~ 1 and E ~ 3 x 10 19 eV, the proton flux can reach the value 1.3 x 
10 24 eV 2 mT 2 sr _1 , which can be compared for example with 2 x 10 24 eV 2 ' mT 2 sr _1 
measured by HiRes. Thus, radiation from cusps may account for observed events at the 
highest energies. The predicted spectrum at E > 8 x 10 19 eV is steeper than the stan- 
dard UHECR spectrum with homogeneous distribution of sources. The accompanying UHE 
gamma radiation is very low, due to large distances between the sources (clusters of galaxies). 

As already mentioned, practically all predicted quantities, such as the diffuse neutrino 
flux (18 ip . the cascade energy density (1821) . the UHE gamma-ray flux from Virgo cluster (1831) . 
the diffuse flux of UHE protons (I8~4"j) and the proton emissivity (I85p . do not depend on the 
basic string parameter r\. There are only two observable quantities that do, the rate of 
neutrino bursts A<b and the neutrino fluence <E>(> E): 

N b ~ lo 2 ^^ y^ 1 ( 86 ) 

IcVlO 
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$(> E) = 1 x 10 



2 fo& AO 10 GW 



) 



1 



-km 



-2 



(87) 



s_ 6 V £ 



[(1 + 3)1/2 _!] 



As ?7 decreases (at a fixed i c ), the rate of neutrino bursts goes up and the number of 
neutrinos detected in a burst, 



We have considered here only " reg ular" , field theory cosmic strings. Recent developments 
in superstring theory suggest [13, IH, [59| that the role of cosmic strings can also be played 
by fundamental (F) strings and by D-branes. Such strings may be superconducting, in 
which case they will also emit bursts of relativistic particles from their cusps. The main 
difference from the case of ordinary strings is that the probability for two intersecting strings 
to reconnect, which is p — 1 for ordinary strings, can be p < 1 and even p <C 1 for F or D- 
strings. A low reconnection probability results in the enhanced density of loops, presumably 
by a factor ~ l/p; the particle production by loops is increased correspondingly. 

UHE neutrinos from superconducting strings may have three important signatures: corre- 
lation with clusters of galaxies, multiple neutrino-induced showers observed simultaneously 
in the field of view of a detector, e.g. JEM-EUSO, and detection of ~ 10 TeV gamma- 
radiation from Virgo, the nearest cluster of galaxies. 
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